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Short dephasing times pose one of the main challenges in realizing a quantum computer. Different 
approaches have been devised to cure this problem for superconducting qubits, a prime example 
being the operation of such devices at optimal working points, so-called "sweet spots." This latter 
approach led to significant improvement of T2 times in Cooper pair box qubits [D. Vion et al., Science 
296, 886 (2002)]. Here, we introduce a new type of superconducting qubit called the "transmon." 
Unlike the charge qubit, the transmon is designed to operate in a regime of significantly increased 
ratio of Josephson energy and charging energy Ej /Ec ■ The transmon benefits from the fact that its 
charge dispersion decreases exponentially with Ej/Ec, while its loss in anharmonicity is described 
by a weak power law. As a result, we predict a drastic reduction in sensitivity to charge noise 
relative to the Cooper pair box and an increase in the qubit-photon coupling, while maintaining 
sufficient anharmonicity for selective qubit control. Our detailed analysis of the full system shows 
that this gain is not compromised by increased noise in other known channels. 

PACS numbers: 03.67.Lx, 74.50.+r, 32.80.-t 



X 



I. INTRODUCTION 

Quantum information processing has emerged as a 
rich, exciting field due to both its potential applications 
in cryptography [l[ and computational speedup 0, 0, 0] 
and its value in designing quantum systems that can be 
used to study fundamental physics in previously inacces- 
sible regimes of parameter space. A promising physical 
paradigm for quantum computers is the superconduct- 
ing Josephson junction qubit 0, 0, 0], which is clas- 
sified into three types according to their relevant de- 
gree of freedom: charge [1, Q , flux [l(| U3 , and phase 
These systems have potentially excellent scalabil- 
ity thanks to well-established fabrication techniques such 
as photo and electron-beam lithography. Unfortunately, 
superconducting qubits currently have coherence times 
which are not yet sufficient for error correction and scal- 
able quantum computation. 

There are several different strategies for enhancing the 
dephasing times in superconducting qubits. One ap- 
proach [131 ] is to improve the properties of junctions and 
materials to eliminate excess sources of 1// noise, whose 
origin remains unclear so far. This is a difficult and costly 
process, but it is likely to benefit a wide range of qubit 
designs when it is successful. A second approach is the 
elimination of linear noise sensitivity by operating qubits 
at optimal working points. So-called "sweet-spot" oper- 
ation has already demonstrated [l4| an increase in de- 
phasing times over previous experiments @ which could 
be as large as three orders of magnitude, and illustrates 
that simple tailoring of quantum circuit design can boost 
qubit performance. In the long run, a combination of 
both strategies will probably be necessary to realize a 
scalable design for superconducting quantum computing. 

In this paper, we follow the second approach and pro- 



pose a new superconducting qubit: a transmission- line 
shunted plasma oscillation qubit, which we call the trans- 
mon. In its design, it is closely related to the Cooper pair 
box (CPB) qubit in Ref. Q. However, the transmon is 
operated at a significantly different ratio of Josephson 
energy to charging energy. This design choice, as we will 
show, should lead to dramatically improved dephasing 
times. 

Two quantities crucial to the operation of a CPB are 
the anharmonicity and the charge dispersion of the en- 
ergy levels. A sufficiently large anharmonicity is needed 
to prevent qubit operations from exciting other transi- 
tions in the system. The charge dispersion describes the 
variation of the energy levels with respect to environ- 
mental offset charge and gate voltage, and determines 
the sensitivity of the CPB to charge noise: the smaller 
the charge dispersion, the less the qubit frequency will 
change in response to gate charge fluctuations. The mag- 
nitudes of charge dispersion and anharmonicity are both 
determined by the ratio of the Josephson energy to the 
charging energy Ej/ Ec- Increasing this ratio reduces the 
(relative) energy level anharmonicity (which limits the 
speed of qubit operations). However, it also decreases the 
overall charge dispersion and thus the sensitivity of the 
box to charge noise. This reduction is important, since 
even with operation at the first-order insensitive sweet 
spot, the Cooper-pair box can be limited by higher-order 
effects of the 1// charge noise [IH, and by the problem 
of quasiparticle poisoning, which can both shift the box 
from its optimal point. 

The transmon exploits a remarkable fact: the charge 
dispersion reduces exponentially in Ej/Ec, while the 
anharmonicity only decreases algebraically with a slow 
power law in Ej/Ec [3. Consequently, by operating 
the transmon at a much larger Ej/Ec ratio than the 
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CPB, one can greatly reduce charge noise sensitivity in 
the qubit while only sacrificing a small amount of anhar- 
monicity. In fact, the charge dispersion can be so strongly 
suppressed that the qubit becomes practically insensitive 
to charge. This eliminates the need for individual electro- 
static gates and tuning to a charge sweet spot, and avoids 
the susceptibility to quasiparticle poisoning, which both 
benefit the scaling to larger numbers of qubits. Amaz- 
ingly, the transmon can at the same time increase the 
strength of electrical coupling between qubits, or between 
a qubit and a transmission line cavity serving as a bus. 

Although the transmon has an Ej /Ec ratio in between 
that of typical charge qubits and typical phase qubits, it 
is important to emphasize that the transmon is very dif- 
ferent from both the CPB and phase qubits, including 
the capacitively shunted phase qubit proposed recently 
by Steffen et al. [13 ■ In the transmon, it is the natural 
anharmonicity of the cosine potential which allows qubit 
operations, whereas in the phase qubit, the Ej/Ec ratio 
is so large that the required anharmonicity can only be 
restored by driving a current I very close to Ic through 
the system, creating a washboard potential, see Refs. [5(, 
@ and Q for recent reviews. The device presented in 
Ref. [ljj operated at an energy ratio of Ej / Ec ~ 2 x 10 4 , 
whereas the transmon will typically involve ratios of the 
order of several tens up to several hundreds and is op- 
erated without the need for any dc connections to the 
rest of the circuit. Thus, the transmon is a new type 
of superconducting qubit that should fix the main weak- 
ness of the CPB by featuring an exponential gain in the 
insensitivity to charge noise. The favorable insensitivity 
of CPBs to other noise sources such as critical current 
and flux noise is maintained (and further improved) in 
the transmon system, rendering it a very promising can- 
didate for the next generation of qubits (see Table |J) . A 
complementary proposal for using a capacitor to modify 
the Ej/Eq ratio in superconducting flux qubits is put 
forward in Ref. [l8j |. 

The outline of the paper is as follows. In Section Hi Al 
we introduce the transmon and its effective quantum 
circuit. The solution of the corresponding Schrodingcr 
equation and an analysis of its asymptotics enable a 
quantitative discussion of the charge dispersion and the 
anharmonicity in Sections III Bl and III Ci respectively. 
Section III El provides additional information about the 
flux degree of freedom in the split transmon, and the 
role of asymmetry in the two Josephson junctions. The 
circuit quantum electrodynamics (circuit-QED) physics 
[l9| of the transmon is investigated in Section IIIII where 
we show that despite the smallness of the charge dis- 
persion, the transmon is expected to reach the strong- 
coupling limit of circuit QED. That is, we show that 
even though the transmon energy levels are insensitive 
to low frequency voltages, transitions between levels can 
strongly be driven by resonant radiation. We discuss in 
detail the modifications of the dispersive limit and the 
Purcell effect due to the increased Ej/Ec ratio. The 
subsequent Sections IIVI and IVl are devoted to the investi- 




FIG. 1: (Color online) (a) Effective circuit diagram of the 
transmon qubit. The two Josephson junctions (with capaci- 
tance and Josephson energy C.j and Ej) are shunted by an 
additional large capacitance Cb, matched by a comparably 
large gate capacitance C g . (b) Simplified schematic of the 
transmon device design (not to scale) , which consists of a tra- 
ditional split Cooper-pair box, shunted by a short (L ~ A/20) 
section of twin-lead transmission line, formed by extending 
the superconducting islands of the qubit. This short section 
of line can be well approximated as a lumped-element capac- 
itor, leading to the increase in the capacitances C g i, C g i and 
C' B and hence in the effective capacitances Cb and C g in the 
circuit. 

gation of noise in the transmon system and its projected 
effect on relaxation (Ti) and dephasing (T2) times. We 
conclude our paper with a summary and a comprehensive 
comparison of the transmon with existing superconduct- 
ing qubits in Section PVTl 

II. FROM THE COOPER PAIR BOX TO THE 
TRANSMON 

A. Model 

In close resemblance to the ordinary CPB (see e.g. 
Ref. Q), the transmon consists of two superconduct- 
ing islands coupled through two Josephson junctions, 
but isolated from the rest of the circuitry. This dc- 
SQUID setup allows for the tuning of the Josephson en- 
ergy Ej = Ej, max |cos(7r<I>/<I>o)| by means of an external 
magnetic flux For simplicity, we initially assume that 
both junctions are identical. (The discussion of the gen- 
eral case including junction asymmetry is postponed un- 
til Section Hi El ) Schematics of the device design and the 
effective quantum circuit for the transmon are depicted 
in Fig. [TJ 

As usual, the effective offset charge n g of the device, 
measured in units of the Cooper pair charge 2e, is con- 
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trolled by a gate electrode capacitively coupled to the 
island such that n g = Q r /2e + C g V g /2e. Here V g and 
C g denote the gate voltage and capacitance, respectively, 
and Q r represents the environment-induced offset charge. 

The crucial modification distinguishing the transmon 
from the CPB is a shunting connection of the two super- 
conductors via a large capacitance Ob, accompanied by 
a similar increase in the gate capacitance C g . As shown 
in Appendix [5] the effective Hamiltonian can be reduced 
to a form identical to that of the CPB system [20| , 



(a) Ej/Ec = 1.0 



H = AE C (n — n g ) - Ej cos (p. 



(2.1) 



It describes the effective circuit of Fig. HJa) in the ab- 
sence of coupling to the transmission line (i.e. disregard- 
ing the resonator mode modeled by L r and C r ), and can 
be obtained from an analysis of the full network of cross 
capacitances as presented in Appendix [AJ The symbols 
n and (p denote the number of Cooper pairs transferred 
between the islands and the gauge-invariant phase dif- 
ference between the superconductors, respectively. By 
means of the additional capacitance Cb, the charging 
energy E c = e 2 /2C s (C* s = Cj + C B + C g ) can be 
made small compared to the Josephson energy. In con- 
trast to the CPB, the transmon is operated in the regime 
Ej > E c . 

The qubit Hamiltonian, Eq. (|2.ip . can be solved ex- 
actly in the phase basis in terms of Mathieu functions, 
see e.g. Refs. @, [Hj]. The eigenenergies are given by 



E m (n g ) — E C a2[ ng+ k(rn,n g )](—Ej/2Ec), 



(2.2) 



where a u (q) denotes Mathieu's characteristic value, and 
k(m,n g ) is a function appropriately sorting the eigen- 
values; see Appendix IB1 for details. Plots for the lowest 
three energy levels Eq, E±, and E2, as a function of the 
effective offset charge n g , are shown in Fig. [2] for several 
values of Ej/Ec- One clearly observes (i) that the level 
anharmonicity depends on Ej/Ec, and (ii) that the to- 
tal charge dispersion decreases very rapidly with Ej/Ec- 
Both factors (i) and (ii) influence the operation of the 
system as a qubit. The charge dispersion immediately 
translates into the sensitivity of the system with respect 
to charge noise. A sufficiently large anharmonicity is re- 
quired for selective control of the transitions, and the 
effective separation of the Hilbert space into the relevant 
qubit part and the rest, 7i = 7i q © Wrest- In the follow- 
ing sections, we systematically investigate these two fac- 
tors and show that there exists an optimal range of the 
ratio Ej/Ec with sufficient anharmonicity and charge 
noise sensitivity drastically reduced when compared to 
the conventional CPB. 



B. The charge dispersion of the transmon 

The sensitivity of a qubit to noise can often be opti- 
mized by operating the system at specific points in pa- 
rameter space. An example for this type of setup is the 
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FIG. 2: (Color online) Eigenenergies E m (first three levels, 
m = 0, 1, 2) of the qubit Hamiltonian (|2.1|l as a function of the 
effective offset charge n g for different ratios Ej /Ec- Energies 
are given in units of the transition energy £01, evaluated at 
the degeneracy point n g — 1/2. The zero point of energy is 
chosen as the bottom of the m = level. The vertical dashed 
lines in (a) mark the charge sweet spots at half- integer n g . 



"sweet spot" exploited in CPBs [21J . In this case, the sen- 
sitivity to charge noise is reduced by biasing the system 
to the charge-degeneracy point n g = 1/2, see Fig. (Ha). 
Since the charge dispersion has no slope there, linear 
noise contributions cannot change the qubit transition 
frequency. With this procedure, the unfavorable sensi- 
tivity of CPBs to charge noise can be improved signifi- 
cantly, potentially raising T2 times from the nanosecond 
to the microsecond range. Unfortunately, the long-time 
stability of CPBs at the sweet spot still suffers from large 
fluctuations which drive the system out of the sweet spot 
and necessitate a resetting of the gate voltage. 

Here, we show that an increase of the ratio Ej/Ec 
leads to an exponential decrease of the charge dispersion 
and thus a qubit transition frequency that is extremely 
stable with respect to charge noise; see Fig.^d). In fact, 
with sufficiently large Ej/Ec, it is possible to perform 
experiments without any feedback mechanism locking the 
system to the charge degeneracy point. In two recent ex- 
periments using transmon qubits, very good charge sta- 
bihty lias been observed in the absence of gate tuning 

Away from the degeneracy point, charge noise yields 
first order corrections to the energy levels of the transmon 
and the sensitivity of the device to fluctuations of n g 
is directly related to the differential charge dispersion 
dEij/dn g , as we will show in detail below. Here = 
Ej — Ei denotes the energy separation between the levels 
i and j. As expected from a tight-binding treatment, 
the dispersion relation E m (n g ) is well approximated by 
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FIG. 3: (Color online) (a) Rotor analogy for the transmon. 
The transmon Hamiltonian can be understood as a charged 
quantum rotor in a constant magnetic field ~ n g . For large 
Ej/Ec, there is a significant "gravitational" pull on the pen- 
dulum and the system typically remains in the vicinity of 
Lp = 0. Only tunneling events between adjacent cosine wells 
(i.e. a full 2-7T rotor movement) will acquire an Aharonov- 
Bohm like phase due to n g . The tunneling probability de- 
creases exponentially with Ej / Ec , explaining the exponential 
decrease of the charge dispersion, (b) Cosine potential (black 
solid line) with corresponding eigenenergies and squared mod- 
uli of the eigenfunctions. 

a cosine in the limit of large Ej /Ec'- 

E m (n g ) ~ E m (n g = 1/4) - ^ cos(2™ 9 ), (2.3) 
where 

e m = E m {n g = 1/2) - E m (n g = 0) (2.4) 

gives the peak to peak value for the charge dispersion of 
the mth energy level. To extract e m , we start from the 
exact expression (|2.2|) for the eigenenergies and study the 
limit of large Josephson energies. The asymptotics of the 
Mathicu characteristic values can be obtained by semi- 
classical (WKB) methods (see e.g. Refs. [U 111 HI). 
The resulting charge dispersion is given by 

(2.5) 

valid for Ej/Ec 3> 1- The crucial point of this result 
is the exp onential decrease of the charge dispersion with 

\J Ej /Ec- 

The physics behind this feature can be understood by 
mapping the transmon system to a charged quantum ro- 
tor, see Fig. [31 We consider a mass m attached to a stiff, 
massless rod of length I, fixed to the coordinate origin by 
a frictionlcss pivot bearing. Using cylindrical coordinates 
(r, ip, z) , the motion of the mass is restricted to a circle in 
the z — plane with the polar angle ip completely speci- 
fying its position. The rotor is subject to a strong homo- 
geneous gravitational field g = ge x in x direction, giving 
rise to a potential energy V = —mgl cos (p. The kinetic 
energy of the rotor can be expressed in terms of its angu- 
lar momentum along the z axis, L z = (rxp)-e 2 = —ih-§-, 



so that the rotor's Hamiltonian reads 
1? 

H I0 t = ^—7^ - mgl cos ip. (2.6) 

Identifying the (integer-valued) number operator for 
Cooper pairs with the angular momentum of the rotor, 
h «-> L z /h, and relating Ej <-> mgl, Ec «-> (h 2 /8ml 2 ), 
one finds that the rotor Hamiltonian is identical to the 
transmon Hamiltonian with n g = 0. 

To capture the case of a nonzero offset charge, we imag- 
ine that the mass also carries an electrical charge q and 
moves in a homogeneous magnetic field with strength Bq 
in z direction. Representing the magnetic field by the 
vector potential A = Bq(— y, x, 0)/2 (symmetric gauge) 
and noting that the vector potential enters the Hamilto- 
nian according to 

p -> p - qA => L z ^ L z + -qB l 2 , (2.7) 

one finds that the offset charge n g can be identified with 
qB l 2 /2h. This establishes a one-to-one mapping be- 
tween the transmon system and the charged quantum 
rotor in a constant magnetic field. We emphasize that for 
the transmon (and CPB) the island charge is well-defined 
so that n has discrete eigenvalues and <p is a compact vari- 
able leading to ip(p) = ^{p + 2tt). In the rotor picture, 
this corresponds to the fact that the eigenvalues of the 
angular momentum L z are discrete and that the "posi- 
tions" ip and p + 2ir arc identical. It is important to note 
that this mapping is different from the tilted washboard 
model used within the context of resistivcly shunted junc- 
tions, see e.g. [27|, and must not be confused with this 
case. 

In the transmon regime, i.e. large Ej/Ec, the dynam- 
ics of the rotor is dominated by the strong gravitational 
field. Accordingly, small oscillation amplitudes around 
p = are favored; see Fig. [31 Perturbation theory for 
small angles immediately leads to an anharmonic oscilla- 
tor with quartic perturbation (Duffing oscillator). [This 
method will be employed in Section III CI to obtain the 
leading-order anharmonicity corrections.] However, the 
charge dispersion e m cannot be captured in such a pertur- 
bative picture. Within the perturbative approach (at any 
finite order) the tp periodicity is lost and the angular vari- 
able becomes noncompact, —00 < ip < 00. Now, in the 
absence of the boundary condition ip(ip + 2ir) = ip(tp) the 
vector potential can be eliminated by a gauge transfor- 
mation. In other words, the effect of the offset charge n g 
only enters through the rare event of a full 27r rotation, 
in which case the system picks up an Aharonov-Bohm 
like phase. This corresponds to "instanton" tunneling 
events through the cosine potential barrier to adjacent 
wells, and explains the WKB-type exponential decrease 
of the charge dispersion. It is interesting to note that 
the nonvanishing charge dispersion is truly a nonpertur- 
bative quantum effect, which can be ascribed to the dis- 
creteness of charge or equivalcntly to the peculiar role of 
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FIG. 4: (Color online) Comparison of numerically exact and 
asymptotic expressions for the charge dispersion and energy 
levels, (a) Charge dispersion |e m | as a function of the ratio 
Ej /Ec for the lowest four levels. The solid curves depict the 
exact results using Mathieu characteristic values, the dashed 
curves represent the asymptotic expansion, Eq. (|2,5[1 , The 
right vertical scale gives the charge dispersion in MHz for 
a transition frequency of 7 GHz. (b) Energy level difference 
Eom = E m — Eq at n g = 1/2 as a function of the Ej/Ec 
ratio. Solid curves show the exact results; dashed lines are 
based on the asymptotic expression (|2,11[) . The vertical scale 
on the right-hand side gives the transition frequencies from 
the ground state to level m in GHz, assuming a charging 
energy of Ec /h = 0.35 GHz. All numerical data are obtained 
for n g = 1/2. 



the vector potential in quantum mechanics leading to the 
Aharonov-Bohm effect. 

The comparison between the exact result for the charge 
dispersion and the asymptotic expansion is depicted in 
Fig. IDJa). The requirements on the largeness of Ej /Ec 
are seen to become stricter for increasing level index. For 
the transmon, we will mainly focus on the lowest two 
levels, for which Eq. (|2.5j) constitutes a very good ap- 
proximation when Ej/Ec > 20. Asymptotically, the 
differential charge dispersion dE i/dn g is dominated by 
the contribution from the first excited level, so that from 
Eqs. (USD and we have 
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7r ei sin(27rn s ) 



(2. 



As a result, the maxima of \dEoi/dn g \ for Ej/Ec = 
20, 50, and 100 are 7.3 xlO" 2 ^, 1.5 xlO" 4 ^, and 8.9 x 
10~ 8 Ec respectively. These values should be contrasted 
with typical values of conventional CPBs operating in the 
limit of Ej/Ec < 1 at the sweet spot n g = 1/2. In that 
case, the charge dispersion relation can be approximated 
by 



Eqi 



[4E c (2n g -l)f 



(2.9) 



At the sweet spot, the system is only sensitive to second- 
order noise, related to the curvature d 2 Eoi/dn 2 of the 



charge dispersion. This is given by (8Ec) 2 /Ej, which 
for Ej/E c = 1 and 0.1 leads to a curvature of QAE C 
and 640-Ec, respectively. A comparison of these numbers 
demonstrates the remarkable robustness of the transmon 
to charge noise. In Section|Vl we will return to this point 
and translate our results into an estimate for the charge- 
noise induced dephasing time. As we will see, the trans- 
ition's dephasing time due to charge noise is exponentially 
increased in the parameter Ej / Ec- Another consequence 
of the drastically reduced charge sensitivit y is that mea- 
surements of the island charge, see e.g. [28], l29j . can- 
not be employed to discriminate the qubit states. Even 
more, Eq. (|2.3[) implies that all higher derivatives of the 
eigenenergies with respect to offset charge become expo- 
nentially small. Thus, more general co ncep ts such as the 
quantum capacitance C q ~ dE 2 /dn 2 [301, which works 
well at the CPB sweet spot, will fail. This simply reflects 
the fact that the inability of charge noise to "measure" 
the qubit state and hence dephase it also means that a 
charge-based qubit readout becomes impossible. Instead, 
we propose a dispersive readout via the cavity which we 
discuss in Section [Til Bl 



C. Anharmonicity of the transmon 

The impressive gain in charge-noise inscnsitivity by in- 
creasing Ej I Ec has to be paid by a loss in anharmonic- 
ity. Sufficient anharmonicity is required to reduce the 
many-level system to a qubit, which ultimately sets a 
lower bound on the duration of control pulses. In the 
following we show that, in contrast to the charge disper- 
sion, the anharmonicity only decreases with a weak power 
law. Therefore, we can find an Ej/Ec range with sig- 
nificantly improved charge-noise inscnsitivity compared 
to the CPB as well as a sufficiently large anharmonicity. 
We define the absolute and relative anharmonicity by 



E\2 — Eqi , 



a/E 01 



(2.10) 



Combining Eqs. (|2.2p and (|2.10[) , one concludes that the 
relative anharmonicity only depends on the effective off- 
set charge and the energy ratio Ej /Ec- In the following, 
we investigate the anharmonicity evaluated at the charge- 
degeneracy point n g = 1/2, so that we can track the full 
crossover from the regular CPB rcgimcjopcrating at the 
sweet spot) to the transmon regime [31| . As shown in 
Fig. [5]Ja), a r then scales as 9(Ej / Ec)^ 1 in the small- 
Ej/Ec limit. For Ej/Ec ~ 9, it changes sign, indicat- 
ing that for larger energy ratios the transition energy Eyi 
becomes smaller than Eq\. The relative anharmonicity 
exhibits a shallow local minimum around Ej/Ec ~ 17.5 
and asymptotically approaches zero for Ej/Ec — > oo. 

The scaling in this limit can be understood in terms of 
perturbation theory in (Ej/Ec)^ 1 <C 1 [13] • Expanding 
the cosine in Eq. (|2.ip around ip = up to fourth order, 
and treating the resulting quartic term in leading order 
perturbation theory, one obtains the following approxi- 
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FIG. 5: (Color online) Relative (a) and absolute (b) anhar- 
monicity at the degeneracy point as a function of the ra- 
tio Ej/Ec- The solid curves show the exact results from 
Eq. (|2.2p . the dashed curves depict the perturbative result 
from Eq. (|2.12[) . (c) Minimum pulse duration r p (blue/gray 
line) and dephasing time due to charge fluctuations T2 (black 
line, Eq. ((5]9}) as a function of E.j/Ec- (Arrows indicate the 
corresponding axis.) Starting in the charge regime, an in- 
crease of the Ej/Ec ratio initially leads to a strong increase 
in the minimum pulse duration, which diverges when the an- 
harmonicity a crosses zero ("anharmonicity barrier"). Above 
the barrier, the operation time scales with a weak power law 
~ (Ej/Ec) 1 ^ 2 - At the same time, the inverse charge dis- 
persion determining the dephasing time due to charge noise 
increases exponentially in (Ej/Ec) 1 ^ 2 - 



achieve, so that the large anharmonicity in CPBs cannot 
actually be exploited fully. Using a typical pulse length 
of 10 ns, we require a minimum anharmonicity of 

|ov in | ~ Op^oi) -1 ~ (10nsx2 7 rxlOGHz)- 1 = 1/200tt. 

(2.13) 

Employing Eq. (|2.12|) . we find that the energy ratio 
should satisfy 20 < Ej/Eq -C 5 x 10 4 , opening up a 
large range with exponentially decreased sensitivity to 
charge noise and yet sufficiently large anharmonicity for 
qubit operations. In other words, the transmon regime 
is reached without paying any serious penalty, and pulse 
generation techniques common for CPB qubits can di- 
rectly be transferred to the transmon qubit. This is fur- 
ther illustrated in Fig. [He) , where the inverse charge dis- 
persion (determining T2 due to charge noise, see Section 
E} and the minimum pulse duration r p are plotted. As 
discussed in detail in Section [V] below, dephasing times 
for the transmon are expected to be significantly larger as 
compared to CPBs. With the projected dephasing times 
of the order of 20 /xs (most likely limited by critical cur- 
rent noise), pulse durations much longer then 10ns could 
be used, making accessible even larger Ej / Ec values and 
greater charge noise insensitivity. 

We emphasize that our considerations regarding t p 
provide a rough and simple order-of-magnitudc estimate 
for the practical Ej /Ec range. A more detailed analysis 
will also have to take into account exact pulse durations 
and shapes [H, Hl| . 



mation for the eigenenergies [see AppendixIClfor details]: 



D. Comparison: transmon and phase qubit 



E m ~ -Ej + ^8E c Ej ( 



where 



Ec 

-t— — I — — (6m 2 + 6m + 3), 
(2.11) 

\/8EcEj/h is also known as the Joscphson 
plasma frequency. A comparison of this approximation 
to the numerically exact result is shown in Fig. Hfb). 
The resulting asymptotic expressions for the absolute and 
relative anharmonicity are 



-Ec 



a r ~ -(8E.j/E c ) 



-1/2 



(2.12) 



depicted in Figs.^a) and (b) [33l ]. 

With these relations and assuming that the transition 
frequency has a value of woi/27r sw 10 GHz typical for 



experiments [341 ], we can estimate the optimal Ej/E, 



c 

range. The resulting absolute anharmonicity is given by 
a = tkooia r . From the frequency spread of a transform- 
limited pulse, we can estimate the corresponding mini- 
mum pulse duration to be t p ~ |a;oiav| 1 . For coherent 
control of the system, the pulse duration must remain 
small compared to T\ and T2. If the total dephasing 
times for the transmon were of the order of a few hun- 
dreds of nanoseconds as in recent experiments on CPBs 
[H]], reasonable pulse durations would be in the range 
of several tens of nanoseconds. It is interesting to note 
that significantly shorter microwave pulses are difficult to 



There exist some remarkable similarities between the 
transmon and typical phase qubits: both operate at 
Ej/Eg ^> 1 (transmon: ~ 10 2 , phase qubit: ~ 10 4 
fl2l . [l7| ) , so that their relevant energy scale is given to a 
good approximation by the plasma oscillation frequency, 
and both are well protected against charge noise. Given 
these similarities and based on the phase-number uncer- 
tainty [37| , one might wonder whether the transmon is in 
fact more closely related to the phase qubit than to the 
CPB. We now show that this is not the case. 

Structurally, the Hamiltonians of the CPB and the 
transmon arc identical, see Eq. (|2.1j) . For both the 
CPB and the transmon the number operator counting 
the charge transferred across the junction is well-defined 
and the phase is compact, i.e. the phase is restricted 
to the interval < <p < 2ir. Increasing the parame- 
ter Ej/Ec smoothly maps the CPB into the transmon. 
By contrast, in the case of a phase qubit there is a dc 
connection between the two sides of the Joscphson junc- 
tion permitting a current (or equivalcntly flux) bias, and 
making the states with phases (p and tp + 2tt physically 
distinct. This topological difference makes it impossible 
to establish a continuous mapping between the transmon 
and the phase qubit via adiabatic changes of Ej /Ec- 

We emphasize that the relationship between transmon 
and Cooper pair box does not imply that the eigenstates 
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FIG. 6: (Color online) Solutions to the qubit Hamiltonian 
(|2.ip in the charge basis. Panel (a) shows plots ol the prob- 
abilities \(n | ip m )| 2 for the presence of n Cooper pairs when 
residing in the transmon eigenstates m — and 1 for three dif- 
ferent Ej/Ec ratios (n g = 1/2). In the limit Ej/Ec > 1, the 
solutions converge to (a discretized version of the) harmonic 
oscillator wavefunctions with increasing width, (b) Fluctua- 
tions of the number of Cooper pairs n as a function of Ej / Ec 
for the first three transmon levels. Solid lines show numeri- 
cal exact results, dashed lines correspond to the asymptotic 
solution (|2,14[) . The Cooper pair number in the ground state 
(first excited state) fluctuates by approximately 1 (2) for an 
Ej/Ec ratio of 100. 



of the transmon are pure charge states. This is illus- 
trated in Fig. [UJa), where the overlap of the transmon 
eigenstates with pure charge states is shown. For in- 
creasing Ej/Ec, the transmon eigenstates spread over 
an increasing number of charge states. However, as de- 
rived in Appendix [Cl and depicted in Fig.[6jb), the charge 
fluctuations only grow slowly as 



1/2 (Ej_\ 1/4 
\8E C J 



(2-14) 



valid in the large Ej/Eq limit. As an example, for 
Ej/Ec = 100 the number of Cooper pairs only fluctuates 
by approximately 1 and 2 in the ground and first excited 
state. In conclusion, the transmon is a CPB operated in 
the Ej/Ec ~> 1 regime with charge Auctions of the order 
of unity. 



fabrication techniques junction parameters vary and typ- 
ically lead to junction asymmetries up to d = fj^qrf^ 1 — 
±10%. In the following, we discuss the effects of this 
asymmetry, which are also known in the context of CPBs; 
see e.g. [14j . 

The case of asymmetric junctions is described by re- 
placing the cosine term in the Hamiltonian (|2.ip by the 
Joscphson Hamiltonian 



Hj = -Ejx COS 01 - E J2 COS (j)2, 



(2.15) 



where 4>i,2 now describe the individual superconducting 
phase differences across the junctions 1 and 2; sec e.g. 
[27| . The usual argument of flux quantization then leads 
to the condition 



2im + 27r$/$o, 



(2.16) 



with integer n, and $0 = h/2e denoting the magnetic 
flux through the SQUID-like ring and the superconduct- 
ing flux quantum, respectively. Defining the effective 
phase difference of the device as ip = {<f>\ +02)/2 and 
Ejs = Eji + Ej2, the Josephson Hamiltonian may be 
rewritten as 

Hj = -Ejy, [cos(7r$/<J>o)cos(^-l-c? sin(7r$/$o) sin0] 



Ejy, cos 




1 + d? tt 



7T$ 
$0 



cos(ip - ip ), 
(2.17) 



where the phase ipo is determined by tan ipo = 
<itan(7r$/$o)- For constant magnetic flux, this phase 
can be eliminated by a shift of variables. As a result, 
our previous results for the symmetric transmon (d = 0) 
translate to the general case by substituting the Joseph- 
son energy by 



Ej -> Ejy cos 




d 2 tan 2 



7T<I> 



(2.18) 



Interestingly, for asymmetric junctions the flux depen- 
dence of ipo may allow for additional qubit control, not 
involving the resonator, by applying ac magnetic fields. 
As with all extra control channels, junction asymmetry 
leads to an additional qubit decay channel from flux fluc- 
tuations, which we will discuss in Section [IVI 



E. Split transmon: The flux degree of freedom and 
junction asymmetry 

In the previous sections, we have ignored the fact 
that the proposed transmon design in fact involves two 
Josephson junctions. Strictly speaking, this is only ap- 
propriate if the two junctions are identical (i.e. they fea- 
ture the same Josephson coupling energy Eji = Ej2). In 
that case the contributions simply add up and our pre- 
vious treatment is valid. However, with current junction 



III. CIRCUIT QED FOR THE TRANSMON 

In close analogy to the situation of the CPB, embed- 
ding the transmon in a superconducting transmission line 
resonator opens up the possibility of control and readout 
of the qubit state - a scenario that has been termed cir- 
cuit QED [l9|, Hi]]. We start from the quantum-circuit 
Hamiltonian for a transmon attached to a superconduct- 
ing transmission line, depicted in Fig. [lja). With the 
Josephson junctions centered in the transmission line, 
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the relevant resonator mode is the I — 2 mode (volt- 
age antinode at the center of the resonator), and it can 
be described by a simple LC oscillator In the real- 
istic limit of large resonator capacitance C r 3> Ce, the 
quantization of the circuit results in the effective Hamil- 
tonian 



H =4Ec(n — n g Y — Ej cos (p + Hu> r aJa 
+ 2peV I ° ms h(a + « t )> 



(3.1) 



see Appendix [A] for the detailed derivation. Here, uj r = 
l/\/L r C r denotes the resonator frequency, and a (cr) 
annihilates (creates) one photon in the transmission line. 
The root mean square voltage of the local oscillator is 
denoted by V® ms = sjfojj r /2C r . The parameter j3 is de- 
fined as the ratio of the gate capacitance and the total 
capacitance, (3 = C g /C^. 

Rewriting the Hamiltonian in the basis of the uncou- 
pled transmon states | i ) , one obtains the generalized 
Jaynes-Cummings Hamiltonian 

H = hy^ujj \j) (j \+f\h) r a ! a + h^2 / g i j \ (a + a f ), 



with coupling energies 

hg lJ = 2PeV? ms {i\n\ ] ) =hg* i . 



(3.2) 



(3.3) 



The general expression (|3.2[) can be significantly sim- 
plified by examining the matrix elements ( i \ h \ j ) , and 
invoking the rotating wave approximation. First, note 
that the asymptotic behavior of the matrix elements 
can be evaluated within the perturbative approach in- 
troduced in Section III CI and detailed in Appendix [C] 
Asymptotically, the number operator assumes the form 
n = -i{Ej/8E c ) 1/4 (b - ft)/V2, so that 

IO- + iH;>|«^(^) 1/4 , (3.4) 



\U + k\n\j)\ g£Zf£^ 



(3.5) 



with |fc| > 1, and b, w denoting the annihilation and 
creation operator for the harmonic oscillator approxi- 
mating the transmon. It is interesting to note that off- 
diagonal matrix elements with an even difference k be- 
tween states fall off exponentially, which can be under- 
stood from the point of view of the parity of the states, 
as well as from the fact that the leading anharmonic per- 
turbation {b + w) A does not mix even and odd states. By 
contrast, matrix elements with odd k > 1 show a slower 
power-law type decay as Ej / Ec — > oo. This is illustrated 
in Fig. [7(a). From Eqs. (|3.4[) and (|3.5|) we conclude that 
nearest-neighbor coupling gi t i±i constitutes the only rel- 
evant coupling in the large Ej /Ec limit. 

Finally, employing the rotating wave approximation to 
eliminate terms describing the simultaneous excitation 
(de-excitation) of both the transmon and the resonator, 



we arrive at the effective generalized Jaynes-Cummings 
Hamiltonian 

H =h^uj \j) {j\+huj r a)a (3.6) 

3 

+ ft5^0 M+ i|i)<i + l|a t + h.c. 



We remark that in contrast to the regular Jaynes- 
Cummings Hamiltonian, this generalized version does 
not allow for an exact analytical solution. The Hamilto- 
nian remains block diagonal. However, due to the pres- 
ence of many transmon levels and nearest neighbor cou- 
pling, the tridiagonal blocks grow in size, and the solution 
of the general case requires numerical methods. The dis- 
persive limit allows for analytical solutions as we show in 
Section HTlBj 



A. The coupling strength of the transmon 

Despite the exponentially decreasing charge dispersion 
for large Ej /Ec, the coupling between cavity and trans- 
mon, expressed by the coupling energies ft<?y, does not 
become small but in fact even increases. This is a central 
message of this paper, and it is crucial for utilizing the 
transmon system as an actual qubit. 

Mathematically, the couplings <?.y are determined by 
a prcfactor containing the capacitance ratio /3, the rms 
voltage of the local oscillator V® ms: and by a matrix el- 
ement of the number operator for Cooper pairs, which 
depends on the energy ratio Ej/Ec- It is interesting 
to note that there is a fundamental upper bound to the 
magnitude of the prefactor. To see this, we rewrite 



(3.7) 



where we have used the relation C r = e r ir /2u> r Z valid 
for a half-wave transmission line resonator, and a = 
e 2 /47reo?ic denotes the fine structure constant. For re- 
alistic values of the effective dielectric constant, this lim- 
its the relative coupling to about 10% of the resonator 
frequency w r . 

While the magnitude of the prefactor is limited, the 
relevant matrix elements between neighboring transmon 
states exhibit an approximate power-law increase as a 
function of Ej/Ec, as depicted in Fig. [7J As a result 
of Eq. (|3.4[) . the couplings g^i+i asymptotically increase 
as (Ej/E c ) 1/4 0. Comparing this with Eq. i[27H]) . we 
find that the increase of gi,i+i can be directly related to 
the increase of the charge number fluctuations. We em- 
phasize that the parameters f3 and Ej /Ec can be tuned 
separately; (3 is essentially determined by the geometry 
of the device, while Ej/Ec can be tuned in situ by the 
external magnetic flux up to a maximum value fixed by 
the device design. 
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FIG. 7: (Color online) (a) Off-diagonal matrix elements of 
the Cooper-pair number operator as a function of the en- 
ergy ratio Ej/Ec (n g = 1/2). Solid curves represent the 
exact result, dashed curves depict the asymptotic behavior, 
Eq. (|3.4[) . The results illustrate that coupling between neigh- 
boring transmon states is the only relevant coupling in the 
limit of large Ej /Ec -(h) Level scheme for the coupled trans- 
mon system. Transmon states are denoted by | i ), i — 1,2, 3 
for the ground, first and second excited state. Photon num- 
bers n p h in the cavity are plotted vertically. The two arrows 
on the | i = 1, n = 2 ) level illustrate the perturbative level re- 
pulsions leading to the dispersive shift. 



This result is quite remarkable. While the sensitivity 
of the transmon spectrum to the dc component of n g de- 
creases exponentially, the ac response to the oscillating 
cavity field increases in a power-law fashion, fn other 
words, the charge dispersion and the magnitude of ac 
coupling are completely disparate. We can illustrate the 
fundamental difference between dc and ac response by 
the following intuitive picture, see Fig. Ojll]. For large 
Ej/Eq, the transmon can be interpreted as a harmonic 
oscillator in the charge basis, with its quadratic poten- 
tial centered at n = n g . Charge noise typically occurs at 
low frequencies so that it can be treated as an adiabatic 
displacement of the oscillator potential (dc response). 
In the general case, this leads to adiabatic changes of 
the qubit frequency and hence dephasing. However, for 
a harmonic oscillator the frequency remains unchanged, 
despite the significant change of the oscillator state un- 
der displacement. Thus, the charge dispersion vanishes, 
and dephasing is exactly eliminated. Remarkably, the 
transmon approaches this ideal point while retaining suf- 
ficient anharmonicity. On the other hand, the question 
of ac response corresponds to driving the oscillator at 
its resonance frequency. Classically, the drive transfers 
energy into the oscillator; quantum-mechanically, it in- 
duces transitions between different oscillator states lead- 
ing to the coupling. This illustrates that strong coupling 
and zero (or exponentially small) charge dispersion are 
in fact not contradictory. This is the central point of the 
transmon: it is highly polarizable and responds strongly 





FIG. 8: (Color online) Mechanical analogy illustrating the 
fundamental difference between dc and ac response, (a) For 
a slow (adiabatic) change in the suspension point of the os- 
cillator, the oscillator mass is displaced but the oscillator fre- 
quency does not change, (b) For an ac drive at resonance, 
energy is forced into or extracted from the system. Quantum 
mechanically, this leads to the coupling between the transmon 
and the cavity field. 



to electric fields at all frequencies. Just as for a har- 
monic oscillator however, the adiabatic response to low 
frequency fields does not lead to changes in the transition 
frequencies. Unlike the harmonic oscillator though, the 
transmon remains moderately anharmonic. 



B. Control and readout: The dispersive limit 

It has been demonstrated that coherent control and 
readout of CPB qubits can be achieved by operating the 
system in the dispersive limit (l9l . [34| . Qubit operations 
are implemented by means of microwave pulses, readout 
corresponds to a measurement of the phase or amplitude 
of the transmitted radiation of a microwave drive field. In 
the following, we demonstrate that these concepts may 
be directly transferred to the transmon: Readout and 
control of the transmon work exactly the same way as 
for the CPB. 

In the dispersive limit, the detunings A.; = uj^i+i — u> r 
between transmon and cavity are large, i.e. in particu- 
lar g> i/|Ao| <C 1, ,goi/|A + a| <C 1. In this case, we 
can eliminate the cavity-qubit interaction to lowest or- 
der by a canonical transformation, sec Appendix [D] for 
details. It is important to note that due to the reduced 
anharmonicity, virtual transitions through excited trans- 
mon states have to be taken into account. Only after this 
can we restrict the transmon Hilbert space to the ground 
state and first excited state. This procedure leads to the 
following effective Hamiltonian 

fko 1 

H cS =-£ k &, + (t^'r + Hx&z) a f a. (3.8) 
Here, the primes signal parameter rcnormalizations: 
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both the qubit transition frequency and the resonance 
frequency of the cavity get renormalized due to the in- 
teraction, ui' r = uj r — X12/2 and uj' 01 = luqi + xoi- (The 
definitions of the partial dispersive shifts Xij will be dis- 
cussed below.) 

The crucial point of Eq. (|3.8p is that the form of 
this Hamiltonian is identical to the dynamical Stark-shift 
Hamiltonian encountered for a CPB coupled to a trans- 
mission line resonator [3, [39[. Remarkably, despite its 
reduced anharmonicity the transmon behaves in a way 
quite similar to a CPB when operated in the dispersive 
regime. This is very convenient as it implies that control 
and readout techniques previously developed for CPBs 
can be transferred to the transmon regime. Specifically, 
the readout proceeds by subjecting the cavity to a mi- 
crowave field close to its resonance frequency. The ac 
Stark effect causes a dispersive shift of the resonator fre- 
quency depending on the qubit state. Consequently, a 
measurement of the phase or amplitude of the transmit- 
ted field is sufficient to infer the state of the qubit. We 
stress that this measurement is very different from mea- 
surements of the quantum capacitance ~ dEf / dn 2 g [3(| , 
which would fail due to the exponentially small charge 
dispersion in the transmon regime. 

The only difference between transmon and CPB re- 
gards the effective dispersive shift \ m Eq. (|3.8p . It is 
given by 



X = Xox - X12/2, 



with 



V ij 



(3.9) 



(3.10) 



and u>ij = ujj — LUi . We emphasize that these expressions 
can be obtained cither through the canonical transforma- 
tion presented in Appendix [Dl or alternatively, through a 
straightforward application of second-order perturbation 
theory. The latter can be understood in terms of pair- 
wise level repulsions between coupled levels as indicated 
in Fig.lftb). 

In contrast to the CPB case, the transmon's disper- 
sive shift consists of two contributions which enter with 
different signs, and which for a pure harmonic oscillator 
exactly cancel each other. The partial cancellation for 
the transmon system with low anharmonicity is compen- 
sated by the increase in the coupling strength g. As a 
result, the magnitude of the transmon's effective disper- 
sive shift is comparable to that of a CPB. 

As we will show now, the contribution of two terms to 
the dispersive shift also leads to interesting new physics 
beyond the usual ac-Stark effect for the two-level case - 
in particular, negative dispersive shifts as well as signifi- 
cantly increased positive shifts, depending on the detun- 
ing of the qubit. Using Eq. (|3.3p . the dispersive frequency 
shifts between neighboring transmon states can be writ- 
ten as 
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FIG. 9: (Color online) (a) Dispersive frequency shift x as a 
function of detuning Ao between cavity and transmon and en- 
ergy ratio Ej/Ec (top panel), for fixed transition frequency 
W01, as well as j3 = 1/2 and n g — 1/2. (b) Cross-sections 
of the 3d plot for Ej/E c = 20 (a), 40 (6), 60 (c), 80 (d). The 
right-hand side y axis shows the dispersive shift in MHz for the 
representative parameters chosen in the text. In the general 
case, the universal axis for h\ (left-hand side) can be trans- 
lated into MHz with the conversion factor -. 1 90 , 1 /?, 7 ^. , ,, • > ■ 

(w l/27r)|(0 I n 1 1 >| 2 

Inset: Level configuration and regimes for the dispersive fre- 
quency shift. 



By combining Eqs. (|3.4p and p. lip , we obtain the fol- 
lowing asymptotic expression for the dispersive frequency 
shift from Eq. valid for Ej/E c > 1: 



1/2 



fix 



Er 



Ec)' 



(3.12) 



The full expression (|3.9p is plotted in Fig.[5J Intriguingly, 
the interplay of the contributions from \oi an d — X12/2 
leads to three distinct regions in the (Ao, Ej /Ec) plane. 
These are separated by the poles of Eq. (|3. 12[) [hA = 0, 
HAq = Ec] where the dispersive limit breaks down. The 
regions are characterized by differing signs for \ [see inset 
of Fig. OUb)] : (i), (ii) for negative detunings, A < 0, and 
positive detunings exceeding the absolute anharmonicity, 
Ao > Ec, the dispersive frequency shift is negative; (iii) 
for small positive detunings < Ao < Ec, which place 
the cavity frequency between the transition frequencies 
ojQi and X becomes positive and obtains compara- 
tively large values. After its special location in parameter 
space, we name this regime the "straddling regime." The 
comparatively large dispersive shift renders this region 
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particularly interesting for qubit operations and readout. 

It is not completely obvious from the outset that the 
straddling regime can be occupied without either violat- 
ing the dispersive constraint goi/ l^ol ><7oi/ l^o + a \ ^ 
1, or losing strong coupling which requires g 3> T = 
max{ft,7}. Here, k and 7 denote the cavity and qubit 
decay rates, respectively. In the following, using realis- 
tic values for the relevant parameters, we demonstrate 
that the straddling regime can indeed be accessed in the 
strong-coupling and dispersive regime, and that it can 
lead to larger overall dispersive shifts \- 

Based on data from a recent experiment [22j], decay 
rates of the order of T/2ir = 2 MHz are achievable. This 
is much lower than the strong coupling g/2ir ~ 100 MHz 
demonstrated in experiments with the qubit operated at 
a transition frequency of uj i/2tt ~ 7 GHz. We argue that 
the charging energy cited in Ref. [1^ can be be further 
lowered to values of the order of Ec/h = 300 MHz. 

First, we consider the case of negative detuning, 
Ao/27r = —150 MHz, i.e. the cavity frequency is higher 
than all transmon transition frequencies. Assuming a 
coupling strength of goi/2ir = 20 MHz, the conditions 
for strong coupling and the dispersive limit are met, 
and we obtain a negative dispersive frequency shift of 
X ~ —1.4 MHz. This should be contrasted with the situ- 
ation of positive detuning Ao/27r =150 MHz in the strad- 
dling regime, i.e. the cavity frequency is located between 
the transition frequencies u)qi and Ui2- Very remarkably, 
this sign change of the detuning leads to an increase of 
the dispersive frequency shift to \ — 3.4 MHz. We have 
also confirmed the validity of the dispersive straddling 
regime by a numerical diagonalization of the full Jayncs- 
Cummings Hamiltonian. 



IV. ESTIMATES FOR THE TRANSMON'S 
RELAXATION TIME (Ti) 

We have argued that a main advantage of the trans- 
mon as compared to other existing solid-state qubits is 
its remarkable insensitivity to charge noise. In this and 
the following section, we investigate in detail the sensitiv- 
ity of the transmon to various noise channels and discuss 
its performance in terms of the projected relaxation and 
dephasing times. 



A. Relaxation by spontaneous emission 

The fact that the transmon qubit couples to the elec- 
tromagnetic field inside the transmission line resonator 
indicates that radiative decay of the transmon is one in- 
evitable relaxation channel. We can estimate the order of 
magnitude of the resulting T\ by a simple semi-classical 
argument, which agrees with the quantum result from 
Fermi's golden rule. The average power emitted into free 
space from an electric dipolc with dipole moment d, os- 



cillating at angular frequency uj is given by 
p 1 ^ 



4iren 3c 3 



(4.1) 



see e.g. [41|. We can obtain an estimate for the trans- 
mon's dipolc moment from the distance a Cooper pair 
travels when tunneling between the two superconducting 
islands. In a recent experiment approaching the trans- 
mon limit [22| . this distance is of the order of L ~ 15 /jm. 
Hence, as an estimate for the dipole moment we obtain 
d = 2eL. As a result, the decay time for the excited 
transmon level due to emission of radiation is given by 



T{ ad = forfoi/P = 



12-Keohc 3 
d?i 



(4.2) 



'01 



For a realistic qubit frequency ujqi/2tt = 8 GHz this leads 
to relaxation times of the order of 0.3 ms and a corre- 
sponding Q factor of 10 7 . 



B. The Purcell Effect 

When a system is placed inside a resonator, its spon- 
taneous emission rate is altered. This effect is known 
as the Purcell effect [42j . It has been observed in mi- 
crowave cavities using Rydberg atoms [43j |. and recently 
also in electrical circuits (23|. For the transmon coupled 
to a transmission line resonator, the same effect will oc- 
cur and each transmon level will experience a different 
change to its spontaneous relaxation rate. In this section 
we compute these rates. 

The simplest way to obtain estimates for the decay 
rates is to apply Fermi's golden rule to the Hamiltonian 
that describes the interaction of the resonator with its 
bath. This Hamiltonian is 



a'bi 



(4.3) 



where bk and b\ are the bath operators for mode fc, and 
Afe determines the coupling strength of the resonator to 
this bath mode. Using Fermi's golden rule and perform- 
ing a continuum limit, the rate for a transition from an 
eigenstate | i ) of the Hamiltonian (|3.6[) to the final state 
l/)is 



lfc,/! A fc / b\,d + a)b k i 



cm; 



(4.4) 

here exemplified for the case of the loss of one photon 
with energy hwk = E- L — Ef to the bath. The reser- 
voir's density of states at energy fttOk is denoted by p(tUk)- 
Defining k = 2Trhp(uj k )\X k \ 2 , this can be written as 



/ K 



«|(/|a|*)r 



(4.5) 
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In the resonant limit of zero detuning, the system's low- 
est excitation consists of an equally weighted superpo- 
sition of the single photon state and the excited qubit 
state. Consequently, we obtain the standard result that 
the resonator's bath will relax this excited state at rate 
«/2. 

By contrast, in the dispersive limit the transmon states 
acquire only a small photonic component (and photons 
only a small qubit component), which can be obtained 
perturbatively from Eq. (|3.6p . 

Kl) = Ki)+ ^ 9j ' j+1 | n -l,j + l) (4.6) 



-^Wln + lj-l). 

UJj—Xj — LJ r 

As a result, the spontaneous emission rate for the excited 
qubit state | 0, 1 ) will acquire the additional Purcell con- 
tribution 



(o,i) = goi 



7* ' = « 



A 2 ' 



(4.7) 



Similarly, in the absence of photons in the cavity the 
higher transmon levels will see a Purcell-induced relax- 
ation rate of 



-,(M+i) _ K g M+i 

Ik ~ I \1' 



(4.8) 



For a cavity with lifetime 1/k = 160 ns (Q = 10 4 ) and 
<7oi /A = 0.1, the Purcell effect leads to a T\ contribution 
of 16 /j,s. 



of robustness with respect to dielectric losses. An ad- 
ditional significant advantage of the transmon geometry 
is that it can be fabricated with single layer processing 
with no deposited dielectric layers that might cause large 
losses. 



D. Relaxation due to quasiparticle tunneling 

The presence of quasiparticles in the system, due to an 
overall odd number of electrons or thermal breaking of 
Cooper pairs, leads to both relaxation and dephasing in 
qubits based on Josephson junctions 0, [45|. Following 
the arguments by Lutchyn at el. [45| we may estimate 
the resulting decohcrence rates for the transmon system. 
With the total number of conduction electrons given by 
N e = nV, n and V denoting the conduction electron 
density and the metal volume respectively, the number 
of quasiparticles may be obtained as 



iV qp = l 



3v/27T 



y/Ak B T A/kgT 

e TP 1 



(4.9) 



Dielectric losses 



valid for temperatures small compared to the supercon- 
ducting gap A. We have assumed there is thermal break- 
ing of Cooper pairs, as well as one unpaired electron, 
which could naturally arise if the finite volume of the 
qubit contains an odd number of electrons. The rate 
of tunneling for one quasiparticle across the junction is 
given by T qp = Sgr/^fi [iH, where S — 1/vV is the 
mean level spacing of the reservoir, V its density of 
states and volume and g F the junction conductance mea- 
sured in units of e 2 /h. Generalizing the expressions from 
j45| to the transmon regime, we obtain for the full relax- 
ation rate due to quasiparticle tunneling: 



Dielectric losses from insulating materials, especially 
amorphous Si02, have recently been marked as a po- 
tentially crucial channel of relaxation in superconduct- 
ing qubits [l3|]. While substrates such as crystalline Si 
and sapphire offer favorably low loss tangents tan S = 
Ree r /Ime r of the order of 10 -6 and 10 -8 at cryogenic 
temperatures, amorphous Si02 has been found to exhibit 
loss tangents as large as 5 x 10 -3 [lj|. Such dielectric 
losses affect the electric fields associated with the qubits 
and cause energy relaxation with a rate proportional to 
tan S. 

These findings clearly call for a cautious choice of ma- 
terials when designing superconducting qubits. In addi- 
tion, we point out that the participation ratio of the dif- 
ferent materials present in the immediate vicinity of the 
actual qubit will play a similarly crucial role. In particu- 
lar, in the transmon design the shunting capacitance Cb 
offers the possibility to accumulate a large percentage of 
the electric fields in a well-controlled spatial region with 
favorable substrates, as opposed to storing them in the 
less well-defined Josephson junction region. As a result, 
we expect good performance of the transmon in terms 



Ti = 1/Ti 



r n 

L qp Jv qp 



^'/ 3 ,n ff ±l/2| e ,n s )| 2 



(4.10) 



The matrix element is the Franck-Condon factor which 
accounts for the shake-up of the transmon collective 
mode due to tunneling of one quasiparticle in the sud- 
den approximation, and in the Ej / Ec 3> 1 limit be- 
comes independent of n g . For our rough estimate, we 
disregard the issue of possible nonequilibrium quasi- 
particle distributions and use typical parameter values 
(V = 300 x 5 x 0.1 ,um 3 , g T = 1, T = 20 mK, aluminum: 
E F = 11.7 eV, n = 18.1 x 10 22 cm~ 3 , v = 3n/2E F ) and 
obtain a relaxation time due to quasiparticle tunneling of 
the order of T\ ~ 1 s. It is important to note that the in- 
creased volume of the transmon does not lead to smaller 
T\ in the quasiparticle channel. In fact, the matrix el- 
ement in Eq. (|4. 10[) leads to an increase of T\ for large 
Ej/Ec because the quantum fluctuations of the coherent 
charge grow larger relative to the charge displacement 
caused by the quasiparticle tunneling. The projected 
temperature dependence of the T\ contribution due to 
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FIG. 10: Number of quasiparticles and contributions to the 
relaxation time T\ due to inelastic quasiparticle tunneling as a 
function of temperature at Ej/Ec = 60. From this estimate, 
tunneling of quasiparticles is not expected to limit the perfor- 
mance of the transmon at cryogenic temperatures. In typical 
dilution refrigerator experiments, (phonon) temperatures are 
of the order of 20 mK, marked in the plot by a vertical bar. 



quasiparticle tunneling is depicted in Fig. 1101 The con- 
clusion we draw from Fig.[TU]is that below 100 mK quasi- 
particles should not lead to significant contributions to 
relaxation in the transmon. This result should be rather 
robust to the actual number of quasiparticles present in 
the limit of T — > 0, since a relevant decrease in T\ is 
expected only when the number of quasiparticles reaches 
several thousands. 



E. Relaxation due to flux coupling 

The coupling of the transmon to an external magnetic 
flux bias allows for an in situ tuning of the Josephson 
coupling energy, but also opens up additional channels 
for energy relaxation: (i) there is an intentional coupling 
between the SQUID loop and the flux bias (allowing for 
the Ej tuning) through a mutual inductance M; (ii) in 
addition, the entire transmon circuit couples to the flux 
bias via a mutual inductance M'; see Fig. [TTJ Here, we 
provide simple order of magnitude estimates of the cor- 
responding relaxation times. For the estimate of relax- 
ation rates due to the mechanism (i) , we assume that the 
overall flux applied to the SQUID ring can be decom- 
posed into the external flux and a small noise term, i.e. 
$ = <E> e + $„ with $„ <C $ e - Then, a Taylor expansion 
of the Josephson Hamiltonian (|2.17[) yields 



Hj^Hj + <S> n A 



(4.11) 



where 



.4 



dH 



9$ 



(4.12) 
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sin I - I cos ip — d cos ( e I sin <p 



As in Section III El E JS 
Josephson energy and d 



7T<I> f 

: Eji + Ej2 denotes the total 
(Eji - E J2 )/Ejyi parameter- 




FIG. 11: (Color online) Model for the estimate of relaxation 
times due to flux coupling, describing (i) flux coupling be- 
tween the transmon's SQUID loop and the external flux bias 
with mutual inductance M, and (ii) flux coupling between 
the transmon circuit and an external flux bias circuit via the 
mutual inductance M' . 



izes the junction asymmetry. Treating the noise pertur- 
batively, one can relate the relaxation rate to the noise 
power spectrum, see e.g. [29j . 



1 



1 

V 2 



|(l|i|0>| 2 M 2 SV,>oi). 



(4.13) 



Here, we have made use of the connection between flux 
noise and current noise determined by the mutual in- 
ductance, S<s> n (uj) = M 2 S[ n {oj). At low temperatures 
ksT <C hujQi the current quantum noise is given by 
Si n (u>) = 2Q(u))hu/R. For a typical junction asymmetry 
of 10% and realistic device parameters (Ej = 20 GHz, 
E c = 0.35 GHz, M = 140$ /A R = 50O) we obtain 
relaxation times ranging between 20 ms and 1 s, where 
the maximum (minimum) T\ is reached for an integer 
(half-integer) number of flux quanta threaded through 
the SQUID loop. 

For the decay channel (ii), we may model the entire 
transmon circuit by a simple LC oscillator with L ~ 
H 2 /4e 2 Ej and C w e 2 /2Ec- Classically, the charge then 
oscillates according to Q(t) = Qa cos cut with oscillator 
frequency oj = 1/ \JLC . Assuming that the energy stored 
in the oscillator is of the order of one energy quantum 
huj, we obtain Qo = \/2Chuj and I(t) = — iosinwi with 
Io = uj\/2Chuj. Through the mutual inductance, this 
oscillating current induces a voltage V[nd(*) = V^smuit 
in the flux bias circuit, where Vq = M'uj 2 V2Chuj. The 
environmental R ~ 50S1 impedance will dissipate the av- 
erage power P = V 2 /2R, which allows for the following 
estimate of the relaxation time: 



hw _ R 
1 P M' 2 w 4 C 



RC 
if 2 



where r\ = M'/L measures the effective coupling strength 
in units of the Josephson inductance. It is crucial to note 
that for the particular case of a SQUID loop and a flux 
bias line exactly centered in the middle of the transmis- 
sion line resonator [cf . Fig. [Hb)] , the mutual inductance 
M 1 identically vanishes for symmetry reasons and relax- 
ation via this channel would not occur. However, when 
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realizing the flux bias line with a coplanar waveguide, it 
is natural to displace the line in order to maximize cou- 
pling to the SQUID loop. The resulting mutual induc- 
tance can be estimated and we obtain values of the order 
of M' = 10<l>o/A. Using realistic numbers for the Joseph- 
son coupling, charging energy, and the environmental 
impedance [Ej = 20 GHz, E c = 0.35 GHz, R = 50ft] 
we obtain relaxation times of the order of 70 ms. There- 
fore, relaxation due to flux coupling is unlikely to limit 
the performance of the transmon system. This can be 
explained by the smallness of the effective coupling 77, 
which for the parameters above is of the order of 10 -6 . 



Coupling to spurious modes, balancing, and 
other noise sources 



It is a common phenomenon that measured relaxation 
(and dephasing) times tend to be shorter than theoret- 
ically predicted, indicating that the microscopic origin 
of decoherence is not yet well understood in all cases 
and that additional, unknown noise channels participate. 
Various candidates for this can be discussed, such as cou- 
pling to spurious resonator modes and pinning and un- 
pinning of vortices etc. In addition, phonon emission 
due to bulk and interface piezoelectricity in solid-state 
qubits has recently been pointed out as another possible 
source for energy relaxation [46| . While bulk piezoelec- 
tricity can be avoided by selecting materials with inver- 
sion symmetry, interface piezoelectricity will always be 
present. With a Rayleigh velocity of vr = 4.92nm/ps 
for surface acoustic waves (SAWs) in Si [47| , one obtains 
acoustic wavelengths of the order of Xr ~ 1 /im at typi- 
cal qubit frequencies. Consequently, there is a significant 
mismatch between the acoustic wavelength and the over- 
all lateral size of the transmon, Xr <C 30 /im. This makes 
an efficient coupling to surface acoustic wave modes less 
likely. In addition, well-defined mesas with step heights 
exceeding the evanescent wave depth will act as mirrors 
for SAWs, which may be further exploited in the de- 
vice design to minimize coupling to SAW modes. An 
accurate estimate of this coupling will require a detailed 
modeling of the transmon's capacitor shapes. Thus, ob- 
taining quantitative estimates for the decoherence rates 
in these additional channels may not always be simple. 
Ultimately, the experiment will determine whether they 
play any role in decohering the transmon qubit. 

We finally emphasize that the design shown schemat- 
ically in Fig. []Jb) may require special care due to the 
possibility of imbalance in the coupling to the transmis- 
sion line resonator. Due to its location between bottom 
line and center pin, the transmon will tend to couple to 
a spurious slotline mode of the resonator in which bot- 
tom and top line are out of phase. It is important to note 
that this issue is eliminated by a design with top/bottom 
symmetry. Alternatively, a tailoring of the relevant ca- 
pacitances in the spirit of a Whcatstonc bridge can com- 
pensate for the problem, thus effectively balancing the 



asymmetry between top and bottom line that is induced 
by the geometric layout. 



V. ESTIMATES FOR THE TRANSMON'S 
DEPHASING TIME (T 2 ) 

Generally, the origin of qubit dephasing can be un- 
derstood in terms of fluctuations of the qubit transition 
frequency due to its coupling to the environment. This 
can be described as noise in the external parameters Aj 
of the Hamiltonian. In the present case of the transmon, 
prominent noise sources are the charge fluctuations as 
well as the fluctuations in the critical current and the 
magnetic flux. We first review the general formalism ap- 
propriate to the treatment of small fluctuations around 
the controlled dc value of the external parameters and 
turn to the investigation of the individual noise channels 
in the subsequent subsections. 

Formally, the qubit Hamiltonian may be written as 



q 2 



=x,y,z 



h u ({\i})a u , 



(5.1) 



and each external parameter can be decomposed into its 
(controlled) dc value and fluctuations around it, Aj = 

X^+SXi. Following Refs. [H and PI, the 

case of weak 

fluctuations can be addressed by a Taylor expansion with 
respect to SXt. To lowest order we have 



huj Q1 



j u—x,y,z 



dX< 



SXj&u + OiSX 2 ), 
(5.2) 

where all derivatives are evaluated at the dc values A' ' . 
The fluctuations <5A^ generally result in two distinct ef- 
fects: (i) For sufficiently low frequencies, the fluctuations 
in the a z component can be treated within the adiabatic 
approximation. In this case, they cause random shifts 
of the transition frequency of the qubit, leading to pure 
dephasing (time scale T2). (ii) Higher frequencies will 
break the adiabatic approximation and induce transitions 
between qubit states (energy relaxation, time scale T\). 
Here, we focus on the dephasing aspect (i), i.e. we con- 
sider low-frequency noise in the a z component. 

It is convenient to define T2 via the law for the decay of 
the off-diagonal density matrix elements, which is given 
by 



iu m t I p -ifg dt'v(t') 



(5.3) 



where v(t) = y~]j ^Rav Assuming Gaussian noise, 

it is straightforward to carry out the noise average, and 
the result can be expressed in terms of the noise power 



/oo 
dT(v(Q)v(T))e- 1 ^ = 
00 



hdX,j 



(5.4) 
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which yields [15j,|48(, 
Poi(i) 



e iW0lt exp 



dui sin 2 (wt/2) 

7T - &v M — / — , >„ 

2-7T (w/2) 2 



(5.5) 



As a consequence, the resulting dephasing critically de- 
pends on the magnitude of the autocorrelation time t c of 
the noise. For correlation times small compared to the 
typical acquisition time, t c <C t, the dephasing follows an 
exponential law 



AnO) 



e M0lt exp 



■\\t\S v (u = 0) 



(5.6) 



The corresponding lincshape is Lorcntzian (homogeneous 
broadening) and the dephasing time T 2 ~ 2/ S v (lo = 0). 
Note that this expression is valid for noise spectra with 
a regular low-frequency behavior. For spectra singular 
at oj = 0, the full equation (|5.5p has to be invoked. In 
particular, this applies to the case of 1// noise, 



2nA 2 



(5.7) 



which has been identified as the typical noise spectrum 
for various noise channels. In experiments, the scaling 
exponent /i in is typically in a range 0.8 < /i < 
1.3. For simplicity, we assume fi = 1 in the following. 
The parameter A determines the overall amplitude of the 
fluctuations and will be specified for the various noise 
sources below; see also Table fl] Denoting the infrared 
and ultraviolet cutoffs by uii and u u , this leads to 



An(i) 



e luJ0lt exp 



a 2 ( dh , 



(5.8) 



valid for W;i < 1 and u u t > 1, see e.g. Refs. |15l.l48ll4£ 



A. Charge noise 

We now relate the differential charge dispersion 
dEoi l ' dn g to the sensitivity of the qubit with respect to 
charge noise. The charge noise observed for transmon de- 
vices in experiments [50( indicates that large fluctuations 
of the offset charge only occur on time scales exceeding 
the typical acquisition time of a single experiment. Addi- 
tionally, small fluctuations may persist within each single 
shot. Both aspects contribute to dephasing and are in- 
vestigated in the following. 

Small fluctuations can be treated in terms of Eq. (|5.8j) 
with the external parameter given by A; — > n g . A rough 
estimate of the resulting dephasing time T 2 is obtained by 
substituting the logarithmic contribution by a constant, 
such that 



T 9 



h 

A 



dE 01 



dn r _ 



Aver] 



(5.9) 



Using Eq. (|2.8|) , we indeed find an exponential increase of 
T 2 for large Ej / E c . Hence, the qubit becomes essentially 
immune to charge noise in this limit. Using realistic pa- 
rameter values [Ej = 30 GHz, E c = 0.35 GHz, A = 10~ 4 
[5lj ] we obtain a dephasing time of the order of T 2 ~ 8 s. 

We now turn to the investigation of slow charge fluc- 
tuations with large amplitudes, which cannot be treated 
within a pcrturbativc scheme. In this case, we explicitly 
write the qubit Hamiltonian as 



hojQi + — cos(27m g + 2Tt5n g (t)) 



(5.10) 



The corresponding decay law of the off-diagonal density 
matrix clement then reads 



Poi(i) 



(5.11) 



e lW0lt ( exp 



. £i_ 

l 2h 



dt' cos[27r(n 9 + Sn g (t'))] 



For variations slow compared to the typical measurement 
time, the effective offset charge will vary for different 
runs but remain constant within each single run, and 
we can substitute n g + Sn g (t) by a single constant. As 
the worst-case scenario, we will assume that the effective 
offset charge randomly switches according to a uniform 
probability distribution on [0, 1]. This results in 



Poi(t) 



/ dn g exp [— it\t cos(27m g )/2fi,] 
Jo 

^'JodexlW- (5.12) 



The envelope of the Bessel function Jo(z) asymptotically 
falls off as \p2pKZ. Thus, using the 1/e threshold as a 
measure for the dephasing time, we obtain T2 ~ e ^fci | ' 
which again increases exponentially with Ej / Ec- For the 
parameter values used above, we find a T2 of the order 
of 0.4 ms. Hence, we find that the increase in the ratio 
Ej I Eq featured by the transmon leads to an exponential 
immunization of the device against charge noise. 

For a full appreciation of these numbers, it is useful to 
contrast our results with the charge dephasing in a regu- 
lar CPB operated at the charge sweet spot. In this case, 
noise is eliminated to linear order and second-order con- 
tributions dominate. Taking these into account, Eq. (|5.9p 
can be generalized and the CPB dephasing time due to 
charge noise may be expressed as [15J 



tt z A 2 d 2 E ai 



Ej 



n g = l/2 



A 2 TT 2 ME 2 ' 



(5.13) 



which for parameter values realized in experiments on 
CPBs 21( leads to dephasing times of the order of T2 ~ 
1 /xs. 



B. Flux noise 

Noise in the externally applied flux translates into fluc- 
tuations of the effective Josephson coupling energy Ej. 
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For simplicity, we consider the symmetric junction case 
d = 0. Since in the experiment flux is used to tune the 
qubit transition frequency via Ej [23| , we consider both 
flux noise at and away from the flux sweet spot. Assum- 
ing that the flux noise is sufficiently small, we can employ 
Eq. (|5.8|) and obtain 



h 
A 



dE 



9$ 



(5.14) 



A 7T 



2EcEjz 



7T$ 7T$ 

sin — tan — 



-1/2 



valid for Ej ^> Eq- In this transmon regime, the device 
is necessarily operated away from points of half-integer 
numbers of flux quanta where Ej = and Ti vanishes. 
For an order of magnitude estimate, we may use a flux 
bias of <& = < I ) o/4 and representative values A = 10 _5 $o 
H, E JS = 30 GHz, and E c = 0.35 GHz. This yields a 
dephasing time of the order of Ti ~ 1 /is. 

It is important to note that Eq. (|5.14[) results in in- 
finite Ti for an integer number of flux quanta. This is 
the flux sweet spot [2l|, where second-order contribu- 
tions [neglected in Eq. (|5.2[) ] dominate. Analogous to 
Eq. (|5.13p . the dephasing time is given by 



T? 



ir 2 A 2 d 2 E Q1 



<9$2 



h$ 2 



$=0 



A 2 7r 4 ^2Ej S E c 



(5.15) 



Using the same representative parameters as above, we 
obtain for the sweet spot a dephasing time of the order 
of T2 ~ 3.6 ms. By comparison, flux dephasing at the 
double sweet spot in the CPB regime yields a dephasing 
time of 



To 



A 2 ir 4 Ejv ' 



(5.16) 



which for realistic parameters gives values of the order of 
1 ms. 



C. Critical current noise 

A second source of fluctuations of the Josephson en- 
ergy consists of noise in the critical current, which is 
believed to be generated by trapping and detrapping of 
charges associated with spatial reconfigurations of ions 
inside the tunneling junction [53[. Such rearrangements 
in the junction directly influence the critical current and 
hence the Josephson coupling energy Ej = I c h/2e. The 
corresponding dephasing time for the transmon is ob- 
tained as 



h 
A 



dE 01 



dl c 



2h 
AE< 



01 



(5.17) 



where A = A/I c denotes the dimensionless fluctuation 
amplitude (independent of Ej). For A = 10~ 6 / c (based 



on [531), E JS = 30 GHz, and E c = 0.35 GHz, this results 
in a dephasing time of the order of T 2 ~ 35 fj,s. For 
comparison, evaluation of the critical current noise in the 
CPB regime yields Ti ~ A ^ , i.e. half the dephasing 
time expected for the transmon. We remark that based 
on our estimates, critical current noise is likely to be the 
limiting dephasing mechanism. 



D. Dephasing due to quasiparticle tunneling 

In principle, the tunneling of quasiparticles not only 
results in relaxation but also in dephasing. For the CPB 
regime, Ej < Ec, Lutchyn et al. demonstrated that the 
dephasing rate may be approximated by ~ r qp iV q p 



[451 ] . where we are using the same definitions as in Sec- 
tion IIVDI A crucial ingredient in this estimate is the 
fact that in the charge limit, the transition frequency of 
the qubit is drastically altered when adding or removing 
a single charge from the island. As a result, in the CPB 
regime complete dephasing is achieved by the tunneling 
of a single quasiparticle. We emphasize that this does not 
hold in the transmon regime, where the charge dispersion 
is exponentially flat so that transition frequency varia- 
tions due to a single charge are minimal. Instead, we ex- 
pect that dephasing due to quasiparticles will mainly be 
induced by tunncling-induced relaxation processes, de- 
scribed in Section ITV PI 



E. Ec noise 

Remarkably, as compared to the Cooper pair box the 
transmon should feature extended relaxation and dephas- 
ing times in all noise channels discussed so far. A possi- 
ble channel for which this does not hold is noise in the 
charging energy, i.e. fluctuations in the effective capac- 
itances of the circuit. So far, there does not seem to 
be any concrete experimental evidence for this type of 
noise. In distinction to the critical current, spatial re- 
configurations of atoms or groups of atoms inside junc- 
tions should only weakly affect the actual capacitances 
(no exponential dependence). Due to the presence of a 
large shunting capacitance in the transmon, Eq noise (if 
existent) could be more important in this system than 
for the CPB. However, presently there is no evidence for 
this to be a limiting factor. 



VI. SUMMARY AND CONCLUSIONS 

In summary, with the transmon we have proposed a 
new type of superconducting qubit: the transmission- 
line shunted plasma oscillation qubit. In terms of the 
value of the energy ratio Ej/Ec ~ 10 2 , the transmon is 
intermediate between a CPB qubit and a current-biased 
phase qubit, but unlike cither of these has no dc connec- 
tions. At the same time, it may be viewed as a natural 
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Noise source 




transmon 


CPB 






Ej/Ec = 85 


Ej/Ec = 1 


dephasing 


1/f amplitude 


T 2 [ns] 


T 2 [ns] 


charge 


A = 1(T 4 - 10" 3 e [51] 


400,000 


1,000* 


flux 


A = 1CT 6 - l(T 5 $o [52, 54] 


3,600,000* 


1,000,000* 


crit. current 


^ = 1(T 7 - 10- 6 / [53] 


35,000 


17,000 



TABLE I: Comparison of dephasing times for the transmon and Cooper pair box qubits. Contributions to T2 are theoretical 
predictions based on Section [V] Entries in bold face mark the dominant noise channel. For the CPB second-order charge noise 
at the sweet spot is most likely limiting the performance of the qubit. In contrast, for the transmon dephasing is suppressed 
to an extent that coherence times should be limited by relaxation (Ti) processes only. All times are given in ns, which is close 
to the clock cycle used in experiments. Typical qubit frequencies are of the order of 1-10 GHz, pulse durations usually range 
between 1 and 10 ns. 

*These values are are evaluated at a sweet spot (i.e. second-order noise). 



improvement of the CPB qubit, which provides its un- 
derlying quantum circuit. 

The realization of the vision of quantum computing de- 
pends crucially on the design of physical systems which 
satisfy the DiVincenzo criteria 55| , in particular the re- 
quirements of sufficiently long coherence time and scala- 
bility. Remarkably, the transmon should offer significant 
improvements with respect to both of these requirements: 
(i) we predict a strong improvement of insensitivity with 
respect to charge noise in comparison to the CPB, see 
Table HJ (ii) the drastically improved charge insensitivity 
should make locking to the charge degeneracy point un- 
necessary, thus simplifying the setup of multi-qubit sys- 
tems significantly. As we have shown, the key to the 
favorable properties of the transmon lies in the combi- 
nation of exponential decrease of the charge dispersion, 
the slow power-law decay of the anharmonicity, and the 
realization of strong coupling to the transmission line res- 
onator. In addition, the dispersive regime of the coupled 
system is described by an ac-Stark shift Hamiltonian in 
complete analogy to the regular CPB, allowing for the 
transfer of control and readout protocols from the CPB 
to the transmon system. 

The effort to reduce the noise susceptibility in solid- 
state qubits based on Josephson junctions has lead to a 
variety of different qubit types. Usually, these designs 
achieve a noise suppression in one particular channel, 
oftentimes accompanied by a tradeoff with respect to 
noise in other channels. Flux qubits [13, [H| operate 
at Ej/Eq ratios similar to those of the transmon, i.e. 
Ej/Ec ~ 10 2 — 10 3 . Accordingly, flux qubits reach an 
insensitivity to charge noise comparable to the transmon. 
However, flux qubits will typically show a significantly 
larger susceptibility to flux noise, especially when oper- 
ated away from the flux sweet spot. Phase qubits [l2j 
trade in a slight increase in critical-current noise sensi- 
tivity for a drastic suppression of charge noise. Recent 
devices using inductive coupling to establish a current 
bias [l7| may also face increased flux sensitivity. 

Remarkably, the transmon achieves its exponential 
insensitivity to 1/f charge noise without incurring a 



penalty in the form of increased sensitivity to cither 
flux or critical-current noise. This advantage can be il- 
lustrated by comparing the transmon to the traditional 
CPB, as shown in TableQ] As discussed above, the trans- 
mon is in fact comparatively less sensitive to flux and 
critical-current noise than the CPB. In fact, even without 
any reduction in the canonical 1// noise intensities, we 
predict that a transmon qubit operated at the flux sweet 
spot should be limited only by the effects of relaxation. 
In conclusion, we are confident that the transmon will be- 
long to a new generation of superconducting qubits with 
significantly improved coherence times and scalability. 
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APPENDIX A: FULL NETWORK ANALYSIS 

For completeness, we describe the modeling of the 
transmon device starting from an analysis of the full 
capacitance network (56[. This network is depicted in 
Fig. [l"27a) . It is based on the capacitances between the 
various conducting islands, see Fig.[T^b). For a minimal 
model, we take into account the two ground planes and 
center pin of the transmission line resonator as well as the 
two islands connected through the Josephson junctions. 
In the actual device, the dc bias is supplied via an addi- 
tional capacitance to the center pin. For simplicity, we 
restrict our network to five islands in Fig. 1121 consider- 



18 



(a) 



HI 



'C35 C34 



C45 



rv 



© 



1-4 



C25 Ci; 



(b) 



(c) 




© CgZ 







FIG. 12: (Color online) (a) Full capacitance network for the 
transmon device, (b) Simplified schematic of the transmon 
device design (not to scale), (c) Reduced network. 



ing only the effective voltage V between bottom ground 
plane and center pin. 

By Thcvcnin's theorem, any single-port linear network 
of impedances and voltage sources can be substituted 
by an equivalent circuit consisting of one voltage source 
V and one impedance. In our particular case it is use- 
ful to retain the original gate voltage source V g in the 
equivalent circuit. This can be accomplished by allow- 
ing for one additional impedance, as shown in Fig. 112( c). 
The two effective capacitances can be interpreted as an 
effective gate capacitance C g and an effective shunting 
capacitance Cb- Together, they adjust for the correct 
voltage seen from the Josephson-j unction port via the 
parameter (3 = V a b/V g = Cg/C^, and the total capaci- 
tance Cs = Cb + C g + Cj between the nodes a and b; 
see Fig. [TSTc). [In the following, we absorb the junction 
capacitance into Cb-] 

The parameters f3 and Cs are extracted from the full 
capacitance network as follows. Each conducting island, 
enumerated by i = 1, . . . , n, is associated with a certain 
charge Qi and a potential <f>i (with respect to infinity). 
These obey the linear relation Qi = J^j Cij<frj- For each 
island, we know either its charge or its potential. Let 
us choose the island enumeration such that for islands 
i < io, the charges Q* are known, whereas for i > io 
the potentials <j>* are known. (Here, the additional star 
signals that the quantity is known.) We thus have the 
following system of linear equations: 

QI = ^2 C ij0j + ^2 C H°i for*<*o (Al) 

j<io j>io 
j<io j>io 

from which we can determine the unknown quantities 



(f>i, . . . , <f>i , Qi +i, . . . , Q n . With the solution, we can im- 
mediately calculate the voltage exhibited at the a—b port 
by subtracting the corresponding island potentials. This 
yields the splitting parameter (3. The total capacitance 
Cs is obtained from the full network by substituting the 
voltage source by a short and calculating the total charg- 
ing energy of the network when applying a voltage Vj 
across the junction. Equating the result with CsVj/2, 
one obtains the parameter Cs- 

The treatment of the transmon embedded in a trans- 
mission line resonator is only slightly more complicated. 
Again, the use of Thevenin's theorem allows for the re- 
duction of the capacitance network to a few effective ca- 
pacitances, see Fig. [Ha). Here, the effect of the resonator 
can be modeled by a local LC oscillator 19!]. Following 
the standard quantization procedure for circuits [20| , we 
obtain 



H = 



4> 2 r | {C B + C g )Ql 
2L r 2CJ 

(Cg + Ci n + C r )Qj 



(A3) 



- Ej cos I 



CgQrQ.J (CsCin + C g Ci n )Q r V g + C g Ci n QjVc 



ci 



ci 



For simplicity, we have absorbed the junction capaci- 
tances into the parallel capacitance Cb, and introduced 
the abbreviation 

= CbC„ + CsCin + C„Ci n + CbC t + C„C r . 



In Eq. (|A3[) the first two terms describe the local oscil- 
lator of the resonator, the two terms in the second line 
capture the qubit's degrees of freedom and the terms in 
line 3 give the coupling between the two of them and the 
coupling to the gate electrode. Taking into account that 



V = V u 



a*) and assuming that C r Cb, Ci n , C g , 



we recover the Hamiltonian (|3.1[) . 



APPENDIX B: MATHIEU SOLUTION FOR THE 
CPB HAMILTONIAN 



We briefly review the solution of the Hamiltonian (|2.ip 
in terms of Mathieu functions, generalizing the results 
from Refs. @, [l6[ to arbitrary values of the effective offset 
charge n g . In the phase basis, the stationary Schrodingcr 
equation is given by 



dip 



Ej cos ip 



ip(cp) = E%jj(ip), 



(Bl) 

where the boundary condition is tf>(ip) = tf)(ip + 2ir). We 
can recast Schrodinger's equation in the standard form of 
Mathieu's equation by introducing the function g(x) = 



-2in a x 



aX tp(2x), so that 



9"(x) 



E E j , . 
— + — C08(2X) 



g(x) = 0. 



(B2) 



19 



The 27r-periodicity of ip(<p) translates into a pseudo- 
periodicity of g(x) with characteristic exponent v = 
— 2(n g — k), where k € Z. Following the notation of 
Meixner and Schafke [57|, Eq. (fB2j) is solved by the 
Floquet-type solution me u (q = — ^^,x). Accordingly 
the eigenenergies E are fixed by Mathicu's characteris- 
tic value, see Eq. (|2.2p . and the wavefunctions can be 
represented as 



c-K.p(in g (p) 



■ me -2[n g -fc(m,n a )] 



Ej_ 

2Er 



(B3) 

The integer numbers k have to be chosen in such a 
way to correctly sort the eigenenergies and eigenstates. 
This implies that k becomes a function of the band index 
m and the effective offset charge n g . Extending Cottet's 
treatment |l6| to cover the full range n g G K, we find 
that this function is given by 



k(m,n g )= y~] [int(2n 9 + £/2)mod2] 



(B4) 



e=±i 



x [int(n 9 ) + £(-l) m {(m + 1) div 2}]. 

Here, int(a;) rounds to the integer closest to x, amodb 
denotes the usual modulo operation, and a div b gives the 
integer quotient of a and b. 



APPENDIX C: PERTURBATION THEORY FOR 
THE LARGE Ej /E c LIMIT 

For completeness, we briefly review the perturbative 
approach employed in Sections III CI and IIIII for large 
Ej / Ec- Starting from the Hamiltonian in the phase ba- 
sis, Eq. (|Bip . one notes that the Josephson energy acts 
as a strong "gravitational force" on the rotor, effectively 
restricting the angle ip to small values around zero. This 
motivates (i) the neglect of the periodic boundary condi- 
tion, and (ii) the expansion of the cosine for small angles. 
Keeping terms up to fourth order, this yields the poten- 
tial energy -Ej + Ejip 2 /2 - Ej(p 4 /2A. 

The Hamiltonian can now be viewed as a harmonic 
oscillator with a quartic perturbation describing the 
leading-order anharmonicity. Due to (i), the "vector po- 
tential" n g can be eliminated by a gauge transformation, 
and the resulting Hamiltonian can be cast in the form of 
a Duffing oscillator 



H = y / 8E c E J (tfb+ 1/2) -Ej 



Ec 
12 



(& + St) 4 , (ci) 



where b, b> denote the regular annihilation and creation 
operators for the harmonic oscillator approximating the 
transmon. The leading-order correction to the eigenen- 
ergies arising from the quartic term is given by 



E 



(i) 



Ec 
12 



(j|(6 + 6 t ) 4 b-) = - T f(6j 2 + 6j + 3). 

(C2) 

Note that in this section, | j ) denotes the pure harmonic 
oscillator state in the absence of any anharmonicity. The 



leading-order correction to the state | j ] 



qp = E C v (»|(S + 6t) 4 |j, 
12 ^ E, - E 4 



(C3) 



i±3 



causes a mixing of | j ) with the states | j ± 4 ) and 
|J±2). 

Noting that for large Ej /Ec the Cooper pair number 
operator can be reexpressed in terms of 



(C4) 



it is straightforward to evaluate the asymptotic expres- 
sions for the matrix elements of n, n 2 involved in the 
charge fluctuations (see Section III D[) and the coupling 
strength (see Section Hll Aj) . 



APPENDIX D: CANONICAL 
TRANSFORMATION 



The elimination of the interaction term in Eq. (|3.2p to 
lowest order in g^+i/Ai (in the dispersive limit) can 
be accomplished by a canonical transformation H' = 
DHD^ . Writing the unitary operator in the form D = 
exp[5 — S^], we find that the appropriate generator for 
the lowest three qubit states is given by 



§ = ^20ta\i + l)(i\ 



Here, the parameters j3i are defined as 

a _ 9i,i+l _ 9i,i+\ 
Pi 



(Dl) 



(D2) 



Employing the Baker-Campbell-Hausdorff relation and 
retaining terms up to order gfj/A 2 , we obtain the trans- 
formed Hamiltonian 

H' = TiWi | i } ( i | + ftw r a)a + hx%,i+i + l 



- hxoitfa | 0) (0 | + ^2 Hxi-i.i - Xi,t+i)tfa \i) (i\ 

i=l 

+ h maa \i + 2)(i\+ h.c. (D3) 



with 



_ ffM+lff"+2[(^+l ~ Uj+2) - (f*H - 

2{U>i+l - Ui - UJ r )(UJ l+ 2 - LU l+1 - U r ) 



and Xij defined in Eq. (|3.10[) . The terms in the last 
line of Eq. I|D3[) describe two-photon transitions, and are 
negligible as compared to the remaining terms due to the 
smallncss of the parameters r\i. 
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